Two-polariton bound states in the Jaynes-Cummings-Hubbard Model 
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We examine the eigenstates of the one-dimensional Jaynes-Cummings-Hubbard (JCH) model in 
the two-excitation subspace. We discover that two-excitation bound states emerge when the ratio 
of vacuum Rabi frequency to the tunneling rate between cavities exceeds a critical value. We 
determine the critical ratio as a function of the quasi-momentum quantum number, and indicate 
that the bound states carry a strong correlation in which the two polaritons appear to be spatially 
confined together. 
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PACS numbers: 42.50.Pq, 03.65.Ge,71.36. 



The investigation of quantum electrodynamics in 
coupled-cavity systems provides insight about the be- 
havior of strongly interacting photons and atoms via a 
variety of interaction schemes [Ij. With the capability 
of tunable coupling and measurement of individual cav- 
ity fields. Coupled-cavity QED could serve as an useful 
tool to address the control of quantum many-body phe- 
nomena d ij. The Jaynes-Cummings-Hubbard (JCH) 
model corresponds to a fundamental configuration ex- 
hibiting the quantum phase transition of light 04111 • 
such a model, single two-level atoms are embedded in 
each cavity and the dipole interaction leads to dynamics 
involving photonic and atomic degrees of freedom, which 
is in contrast to the widely studied Bose-Hubbard model. 
Recently, the phase diagrams and excitation spectra of 
the JCH model have been discussed in literature [l^ - fTsj . 

In this paper we examine the eigenstates of the JCH 
Hamiltonian in the two-excitation subspace. Our focus 
is placed on the existence of bound states as well as 
their features. It is interesting to note that two repulsive 
bosonic atoms can form a bound pair in an optical lattice 
[l6l - [T8| . As we shall indicate below, the JCH Hamilto- 
nian also supports two-polariton bound states when the 
photon-atom interaction is sufficiently strong. In partic- 
ular, the two polaritons associated with the bound states 
exhibit a strong correlation such that they stay close to 
together in position space. 

To begin with, we first specify the system configura- 
tion and the model Hamiltonian. The one-dimensional 
JCH model consists of a chain of A^-coupled single-mode 
cavities and each cavity contains a two-level atom. The 
cavity mode frequency is uJc and the atomic transition 
frequency is Wa- The Hamitonian of the model in the 
frame rotating with the atomic frequency Wa is given by 
(S=l): 
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the annihilation and creation operators for photons of 
the n-th cavity. The cr^ are Pauli ladder operators for 
the atom in the n-th cavity. The J describes the photon 
tunneling rate between neighboring cavities, and g is the 
vacuum Rabi frequency which characterizes the photon- 
atom interaction strength. In this paper we assume the 
periodic boundary condition such that the cavity labeled 
by n = A'^ -f 1 corresponds to the cavity n — \. 

Defining the atomic and photonic excitation number 

operators as Na = Yll^i^t^n and iV^ = Y.'^=i°'n°'^^ 
it is easy to check that the excitation number (or the 

polariton number) is a conserved quantity, i.e., [H, Na + 
Nj] =0. In this paper we will focus on states with 
two excitations only, i.e., {Na + Nj)\^jj) = 2\^). 

In order to exploit the translational invariance, we de- 
fine the following operators via discrete Fourier trans- 
form: 
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where k — 0, 1, 2, — 1 is related to the (discrete) 
quasi-momentum of the excitation. The commutation 
relations for the operators are given by: 



Sfc,S 



= 5. 
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where tr^ is the Pauli-Z matrix for the n-th atom. The 
Hamiltonian in term of operators defined in ([2]) and ([3]) 
is partially diagonalized: 
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where fife = A -I- 2 J cos are normal mode frequencies 
where A = — is the detuning. The a„ and aj^ are (shifted by Wq) of the field in the coupled-cavity system 
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in the absence of atoms. For later purpose, the ground 
state of the system is denoted by |<i>o) which is the state 
with vacuum cavity fields and all atoms being in their 
ground levels. 

The two-polariton subspace can be spanned by the 
kets: \kj)F = 6l6]|*o), \kj)A = 441*°)' ^^'^ 
\k)F\j)A = ^141*^0^ with k and j ranging from to 
iV — 1. The subscripts F and A are used to denote 
the field and atomic excitations respectively. For conve- 
nience, \kj)F and \kj)A are not normalized. Note that 
\k j)A are generally not orthogonal to each other because 
{k' f\kj)A = Skk'Sjj' + SkjSkk'Skj' - jfSk+j,k'+j' ■ Fur- 
thermore, it can be shown that 



his]\kj)F = 5ij\k)F\i)A + 5ik\j)F\k)A 



(7) 
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{b\si 4- bis\)\k) f\3) A = Sij\kj)F + Sik\kj)A 
b\si\kj)A = 5ki\k)F\i)A + 5ji\j)F\k)A 

-^\l)F\[k + ]-l])A (9) 

with [x] = X (mod N). Eqs. ([7][9l) imply that when 
H operates on |fcj)_F, \k)F\j)A or \ki)A, the quantum 
number P = k~\-j (mod N) remains unchanged. We may 
call this as a conservation of quasi-momentum which is 
the key to construct eigenvectors. 

A general two-polariton eigenvector with a given quasi- 
momentum quantum number P is given by: 

\^p) ^ {akj\kj) F + Pkj\k) f\j) A 

+f3'k^\j}F\k)A+lk,\kj}A) (10) 

where Sp denotes the set of (fc, j) satisfying k + j = P 
(mod N) and j > k. To avoid double counting, we set 
l^'kk ~ 0- Next by the Schrodinger equation, assuming A 
is the eigenvalue, H\'^p) = A|4'p), we have, for j > k, 

Xakj = {nk + n^) akj + g {Pkj + f3'kj) (11) 

X/3kj = gakj + ^kPkj + gikj - ^ X! ^^'3' (12) 

Sp 

>^Pkj = go^k] + ^jP'kj + gikj - ^ XI (1^) 



^Ikj = 9{hj+P'kj) 
and for j — k, 

Xakk = "^^kockk + gPkk 



(14) 
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A/3fcfc = 2gafcfc -I- VLkl^kk + '^.g-fkk - '-^'^Ik'j' (16) 



A7fcfc — gPkk- 
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It is worth noting that the case of even N and odd P is 
simpler because this excludes the possibility k — j. 

We may cast Eqs. ([TTI -17) in the matrix form and solve 
the eigensystem directly by standard numerical packages. 



Since we are interested in systems with a large TV ^ 1, 
the corresponding eigenvectors are expected to be insen- 
sitive to parity of N and P and this has been verified in 
our numerical calculations. To facilitate our discussions, 
we shall focus on the case of even N and odd P so that 
Eq. (15-17) are not needed, and Eqs. (11-14) become: 



HpV = AV 



where Hp is a 2N x 2N matrix: 
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where {ki,ji) G Sp, and h^.j-^ and wo are 4x4 subma- 
trices: 
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In this way the eigenvector V takes the form: 



V= (vfciji,Vfc,,j,,--- ,Vfc 



N/2 ^3 N/2 J 
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with Vfcj = (ctkj, Pkj, fi'kjTlkj'^ ■ We remark that Hp 
is not symmetric and this is due to the fact that non- 
orthogonal atomic basis vectors \kj)A have been used to 
express the eigenvectors. 

In Fig. 1, we illustrate the eigenvalues of as a func- 
tion of P at resonance (A = 0) for various photon-atom 
interaction strengths. We find that if g/J is sufficiently 
large (Fig. Ic and Id), then there exist two branches of 
discrete eigenvalues (marked in red) . The discrete eigen- 
values correspond to bound states in the large N limit 
and they are isolated from the quasi-continuous bands 
of eigenvalues. By lowering the photon-atom interaction 
(Fig. lb), the bands and the discrete eigenvalues become 
closer, and eventually the discrete eigenvalues disappear 
in the bands if g is too small (Fig. la). We will estimate 
the critical value of g for the occurrence of bound states 
near the end of this paper. 

Let |Af,) be the normalized eigenvector with the dis- 
crete eigenvalue A;,. Although the notion bound states is 
for infinite N systems, the confined nature of |Af,) starts 
to emerge in real space at finite large N. This is illus- 
trated in Fig. 2 in which the probability distribution 
of excitations in real space are plotted. Specifically, we 
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FIG. 1; (Color online) Eigenvalues of the H ior N — 50 and 
odd P and A = for (a) g/J = 0.1, (b) g/J = 1.7, (c) 
g/J = 2 and (d) g/J — 5. The eigenvalues are joint by lines 
for ease of visualization. 



calculate the joint probabilities defined by, 
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FIG. 2; (Color online) Joint probabilities distribution asso- 
ciated with bound states for the system with A = 0, P = 1 
and A'^ = 50 at various coupling strengths. The right column 
corresponds to g slightly above the critical value. 



Here m and n are indices for the cavity positions, and 
hence [Pnm) is the joint probability of having pho- 
tonic (atomic) excitations in the n-th and m-th cavity. 
Similarly, is the joint probability of a single photon 
in n-th cavity and an excited atom in the m-th cavity. 
For a given value of P, all the three joint probabilities 
depend on n — m only. 

In Fig. 2(a-c), we see that the joint probabilities are 
strongly localized around m = n when the interaction 
strength g is sufficiently away from the critical value 
(.9c ~ \/3J). Note that p^^^ vanishes at n = m because 
(T+^ = 0, but p^^ peaks at n — m ± 1. The localized 
feature means that the two polaritons are spatially con- 
fined together and behave like a composite object. As 
the photon-atom interaction g decreases, the binding be- 
comes weaker. This is shown in the wider distributions 
in Fig. 2(d-e) in which g/J is slightly above the critical 
ratio. We remark that although Fig. 2 is plotted with 
P = 1 as an illustration, our numerical results indicate 
that the localized feature of bound states appear for the 
whole range of P. For example, in the case = 50 and 
g/J = 5, by adding up the same site probabilities and 
the nearest neighbor probabilities, the sum is over 98.5% 
for < P < A^- 1. 

Although an exact analytic solution of the bound 
states is not available, some analytic insights may be ob- 
tained in certain limiting cases. From Eqs. (|llj|14[) . non- 
zero eigenvalues A are the roots of an algebraic equation 



G(A) = 0, where G is defined by. 
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Cfc,(A) = (2A - rife - %) (A - - ^j) 
Dkj{X)^g' (2A-r!fe-%)' 

-A (A - flk) (A - n^) {X^flk- Qj) .(28) 

In the strong coupling limit with g ^ J, we find that the 
eigenvalues A;, for bound states are approximately given 
by, 
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This is obtained by noting that at J = (i.e., when the 
cavities are decoupled), the iV^g are dressed energies 
known in the resonance Jaynes-Cummings model. Treat- 
ing this as a zeroth order approximation, we can make a 
Taylor expansion of Eq. ([25]) in a power series J/g. Eq. 
(|29p is then obtained as an approximation solution of Ab 
up to the first order in J/g. 

In the ^ 1 limit, Eq. (12^ can also be used to 
estimate the critical value of g below which bound states 
cease to exist. This is done by observing that the function 
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FIG. 3: (Color online) An illustration of the function G for the 
system with A = 0, iV = 50, P = 15 and g/J = 2. Non-zero 
eigenvalues of the Hamiltonian are roots of G(A) = 0. The red 
dots inside the gaps are roots corresponding to bound states. 
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FIG. 4: The critical ratio Qc/ J as a function of P for A = 
systems. 

G(A) form bands, and Af,'s exist in the gaps between the 
bands (Fig. 3). As g decreases the gaps between bands 
become narrower. At a critical g — gc, the bands start 
overlap and no bound states are supported for g < gc- 



Therefore gc can be estimated by determining when the 
band edges overlap. In the I limit, the bands are 

filled up by those A's that give zero denominators in (|26|) 
[i.e., G{X) — ±oo] as shown in Fig. 3, and hence the band 
edges can be estimated. Base on this estimation scheme, 
the critical value gc in the A^ — > oo limit is plotted as a 
function of P (Fig. 4). Note that gc is a maximum at 
P = and P = TV - 1 and it equals %/3 J. 

To conclude, we have addressed the spectrum as well 
as the conditions for the existence of two-polariton bound 
states of the JCH model. These bound states are compos- 
ite objects characterized by the quasi-momentum quan- 
tum number P and they process a strong spatial corre- 
lation. In view of the previously studied bound atom 
pairs in Bose-Hubbard model [iGl - flsj , our work is a gen- 
eralization to bound polariton pairs in a coupled cavity 
QED system. Although our analysis have been confined 
to the zero detuning case, bound states are also found for 
detuned systems according to our numerical solutions of 
Eq. (18). The effect of detuning would control the na- 
ture of polaritons such that the two excitations can be 
mainly photonic (atomic) at positive (negative) large A. 
We also remark that the dynamics of a single pola riton 
in the JCH model has been discussed recently [1^, and 
we expect the two-polariton problem may have richer dy- 
namical features because bound states would enable in- 
teresting correlated two-polariton transport. We hope to 
address this issue in the future. 
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